Complementary cosmological tests of RSII brane models 
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In this paper we explore observational bounds on flat and non-flat cosmological models in Type 
II Randall-Sundrum (RSII) branes. In a first analysis, we consider current measurements of the 
expansion rate H (z) (with two priors on the local Hubble parameter) and 288 Type la supernovae 
from the Sloan Digital Sky Survey (within the framework of the mlcs2k2 light-curve fitting method). 
We find that the joint analysis involving these data is an interesting tool to impose limits on the 
brane tension density parameter (fix) and that the spatial curvature has a negligible influence on 
fix estimates. In order to obtain stronger bounds for the contribution of the fix we also add in 
our analysis the baryon oscillation peak (BAO) and cosmic microwave background radiation (CMB) 
observations by using the so-called CMB/BAO ratio. From this analysis we find that the fix 
contribution is less than 4.10 -5 (If). 

PACS numbers: 98.80.-k, ll.25.-w 

I. INTRODUCTION 

Since the formulation of Kaluza-Klein theory in the beginning of last century, concerning the unification of gravity 
and electromagnetism in a five-dimensional space, the study of extra dimensions have attracted much attention 
in Physics. Nowadays there are great interest and much speculation about the existence of extra dimensions and 
their influence over our four-dimensional world, motivated mainly by the braneworld scenarios, in which our four- 
dimensional spacetime is viewed as a submanifold isometrically embedded in a space of higher dimensions [J 0, y, l3, 

a a. 

A relevant aspect of the braneworld scenarios is that matter and all fields are confined to the brane(below a certain 
energy level which is expected to be of Tev order) and only gravity can propagate into the bulk [H 0, S Hi • Because 
of this characteristic, extra dimensions might have large scale as compared to the tiny Planck scale, assumed by the 
Kaluza-Klein theory, without any conflict with the observations so far. In particular, the RSII model admits the 
existence of a hidden extra dimension which has an infinity length |^|. If the hidden extra dimensions are large, 
in principle, they can be detected more easily by exploring the signature they might leave in our four-dimensional 
space. This could be done by trying to detect deviations of four-dimensional laws in domains where gravity was not 
empirically tested yet. Pursuing this goal, some research teams have investigated the validity of the inverse square law 
of gravity at sub- millimeter scales recently 0, S 3 ■ According to them, there is no evidence of extra dimensions down 
to a length scale of about 50/zm Q . In cosmology, there are also many studies about the impact of extra dimensions 
in the cosmic evolution considering distinct braneworld models [y, [ljj, [ll|. Here we are interested in discuss some 
cosmological tests of RSII model. 

The RSII model is distinguished from other braneworld models because it is the first one which, surpassing the 
paradigm of compact extra dimensions, shows the possibility that we live embedded in a space with an extra dimension 
of infinite size. Despite the fact that gravity can have access to the infinite extra dimension, a negative cosmological 
constant A( 5 ) defined in the bulk ensures the existence of a massless graviton mode, which is confined to the brane, and 
the suppression of the massive modes. These two characteristics together, on their turn, guarantee that the Newtonian 
law is recovered at distances much larger than the curvature scale £ = ^/— 6/A( 5 ) established by the five-dimensional 
cosmological constant [J]. 

The laboratory tests of the inverse square law, that have found no trace of extra dimension, as we already mentioned, 
then suggest that the curvature scale I could not be greater than 10 _4 m [9|. With such a small value, it seems that 
the extra dimension will play no significant role in the recent phase of evolution of the Universe, which, according to 
Type la supernovae (SNe la) data and based on the LCDM model, is undergoing an accelerated expansion. Needless 
to say, however, that any theory should be tested in different scales and by means of all kind of physical phenomena 
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as much as possible. Hence, in this paper, we intended to use a recent SNe la data combined with measurements of 
H(z) as an independent and new test for the RSII model in the cosmological scale. 

As is well known, in the cosmological context, the differences between the RSII model and the standard cosmological 
model can be synthesized as modifications of the Fricdmann equation. Admitting that the brane is a Roberton- Walker 
three-dimensional space, then the effective field equation of gravity induced on the brane yields the modified Fricdmann 
equation |f0| : 

£)'- n * + *-(£) +n "(T) ,+n »(£)' +n -(?) 4 o 

where a is the scale factor, H is the Hubble parameter, p is the density of the dust matter confined to the brane. 
The number zero as a subscript indicates the current value of the quantities. The parameters f^A, H m and Q K are the 
density parameters of the four-dimensional cosmological constant A( 4 \, of the matter and of the spatial curvature, 
respectively. They have the usual definition J7a = A^/3H§, Q m = 8ttGpq/3Hq and fl K = —k/oqHq. The new terms 
are: the density parameter of the positive brane tension (A) 

48nGp 2 

and the density parameter of the so-called dark radiation Q u — m/a^H 2 , where m is another free constant of the 
model. 

It is important to emphasize that all terms of Eq. (TTJ are quantities that inhabit the brane, except the dark radiation, 
which is essentially a fivc-dimcnsional quantity and, therefore, represents a direct influence of the bulk geometry over 
the brane [6J, |lOj. The physical origin of the dark radiation can be understood by examining this picture from 
the perspective of the embedding formalism. It can be shown that a homogeneous and isotropic three-dimensional 
space, i.e., a Robertson- Walker 3-space, can be isometrically embedded into the AdSs-Shwarzschild spacetime p^.fl3|. 
which is a black hole solution in five dimensions with a negative cosmological constant. From this point of view, the 
expansion of the observable Universe is a consequence of the motion of the brane in the static AdS5-Scwharzschild 
space. Moreover, it can also be shown that the trajectory of the brane is dictated by equation (fTJ), where the scale 
factor localizes the brane in the bulk and m is the mass of the bulk black hole [12|, [l3| . 

Other remarkable difference of the equation ((T|) , in comparison with the usual Fricdmann equation, is that it depends 
on the energy density squared. The coefficient of this quadratic term is the parameter density Q\ which is inversely 
proportional to the brane tension. So, the higher the brane tension the lower the effects of the extra dimension in our 
world. This happens because, in the RSII model, the brane tension is connected to the five dimensional cosmological 
constant A( 5 ). A high tension corresponds to a high cosmological constant and, therefore, a strong suppression of the 
massive graviton modes. Of course measurements of the brane tension A by estimates of fl\ would give us important 
information about the influence of extra dimension in the cosmic evolution. 

At this point, it is interes ting to highlight that p 2 arises too in cosmological models based on other theories as the 
Loop Quantum Cosmology [14j and non-Ricmannian theories |15| . Therefore, although the RSII brane model is the 
main motivation of our discussion, the constraints on the density parameter fl\ we find here could also be used to 
put bounds on the parameters of those theories. 

The modified Fricdmann equation depends on the five density parameters that should be determined by empirical 
data. The first confrontation of RSII model with SNe la data was done in Ref. [Tg|. The authors found, by using 
the Pcrlmuttcr samples 17], that SNe la data provide weak constraints on the brane tension. The simplest case with 
rn = and with only dust matter on the brane yields 0^ = —0.9, fl\ = 0.04, tt m = 0.59 and Ha = 1-27 as the best 
fit according to the \ 2 method. Considering these numbers as an unrealistic result, the prior fi m = 0.3 was added in 
the flat model. In this case, H\ < 0.037 at 2a level, by using the Perlmutter sample A jig . 

In Ref. [18( , another data set of SNe la was employed in order to constrain the parameters of the RSII model. The 
best fit, in the flat model with the prior O u = ± 0.1, is H\ = 0.026, Cl m = 0.15 and fl\ = 0.80, by using data of the 
Riess gold samples [ljj. Once again, the constraints were not so strong. For instance, Q m <G [0, 0.62] at 2a level. The 
curved model was also investigated, however, in this case, it was necessary to impose some priors on density matter 
(H m ~ 0.3) in order to get reasonable values for the density parameters. 

The contribution of the p 2 in cosmology was again considered in Ref. [2fJ , this time by using a data set of SNe la 
containing samples from [19|, ]2M, 122] ■ In that paper, the authors also considered other tests like CMB and BAO in 
order to constrain the density parameter of the RSII model. However, taking into account only the SNe la data, they 
found similar results. Indeed, without imposing any additional prior, the best fit obtained (Hk = 0.34, O m = 0.0 , 
£l\ = 0.044 and Ha = 0.646) was again not very reasonable. 

Therefore, based on these previous results, it seems that, as SNe la data do not constrain strongly the density 
parameter of the matter, we do not find realistic results without assuming some priors for the matter density. In this 



paper, we want to show that this can be avoid if we consider a joint analysis with SNc la data and H (z) measurements. 
We explore flat and non-flat cosmologies. As we shall see, with only 21 data points, H(z) test gives supplementary 
information which naturally put realistic bounds on the parameters without manipulating priors to the matter density. 

For the sake of simplicity, we are going to admit hereafter that m = 0. Of course, this case corresponds to the 
simplest model, in which there is no black hole in the original bulk and, therefore, no dark radiation in the brane. 
Nevertheless, we have to mention that there are some papers that have focused their attention on the analysis of dark 
radiation in RSII branes also by using SNe la data, see for example [23J. 

In order to obtain stronger bounds for the parameters of the RSII model, we have also considered the CMB and B AO 
cosmological tests. However, differently from [2(j, our analysis is based on the CMB/BAO ratio which is expected to 
be weakly cosmological model dependent according to [24[ . 

The paper is organized as follows. In section II, we give a short description of the observational data we have used. 
The corresponding constraints on the cosmological parameters are investigated in section III. The article is ended 
with a summary of the the main results in the conclusion section. 

II. DATA SETS 

A. SNe la Sample 

In our analyses we have considered the SNe la data sample from the Sloan digital sky survey [24|, (25| ■ This sample 
is a combination of 103 SNe la with redshifts 0.04 < z < 0.42, discovered during the first season of the Sloan Digital 
Sky Survey-II (SDSS-II), with a comprehensive and consistent reanalysis of other datasets [ljj, |2& [27[ , resulting in a 
combined sample of 288 SNe la. These 103 data filled the redshift desert between low- and high-redshift SNc la of 
the previous SNe la surveys. Kessler et al. (2009) used two light-curve fitters to obtain the SNe la distance moduli, 
namely, MLCS2K2 [H| and SALT2 [H. In this paper we use the SNe la distance moduli obtained via MLCS2K2 
calibration since this method does not have cosmological model dependence [3fJ|. Furthermore, as is largely employed 
in literatura, all the results in our analysis (see next section) from SNe la data are derived by marginalizing the 
likelihood function over the nuisance parameter [13 . |31| . 



B. H(z) measurements 

In recent years H (z) measurements of the Hubble parameter have been used to constrain cosmological parameters 
[22, [33, [34|, [35|, [3a, [33. The idea underlying this approach is based on the measurement of the differential age 
evolution of the massive and passively evolving early- type galaxies as a function of redshift, which provides a direct 
estimate of the Hubble parameter H(z) — —1/(1 + z)dz/dt « —1/(1 + z)Az/At. It is important to stress that direct 
measurements of H(z) at different redshifts is also possible through measurements of the radial component of baryonic 
acoustic oscillations [38| (for recent H(z) reviews see [331]). 

In this paper, we have used the 21 H(z) measurements [34|, [39, Ho] in the redshift range < z < 1.75. The H 
influence on our results is explored by considering two priors in analysis: 68 ±2.8 km/s/Mpc, from a median statistics 
analysis of 553 measurements of Hq [35|, and 73.8 ± 2.4 km/s/Mpc, the most recent one based on HST measurements 
[41j . Following J37[ we have assumed that the distribution of Hq is a Gaussian with one standard deviation width 
a h q an d mean Hq, so that we can build the posterior likelihood function Cn{p) that depends only on the parameters 
p by integrating the product of exp(— Xh/^) an< ^ the -^0 prior likelihood function exp[— (Hq — Ho) 2 /(2a H )], 

00 

Ch (p) = 1 I e -X 2 H(^,p)/2 e -(ff„-S„)V(2<) rfiJo , (2) 

In this way, to obtain our results from H{z) measurements we maximize the likelihood Ch(p), or cquivalently minimize 
Xh(p) = — 21n£fl-(p), with respect to the parameters p to find the best-fit parameter values to flat and non-flat 
universes (see next section). 

C. Ratio CMB/BAO 

As it is largelly known, the CMB constrains from the so-called shift parameter R [42( and the BAO measurement 
A [43J have been commonly used to constrain non-standard models, but this approach is not completely correct since 
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FIG. 1: Confidence contours on the plane Q\ — Qm to flat Universe. Fig. la) Contours are drawn to la and 2a. The black 
dots corresponds to bounds from SNe la. Through this paper the H(z) priors 68 ± 2.8 and 73.8 ± 2.4 correspond to blue solid 
and red dash-dot curves, respectively. Fig. lb) shows the contours to la and la from our joint analysis. Fig. lc) displays the 
likelihoods to Q\. 

these quantities are derived by using parameters close to standard wCDM [24|, |42, Hi • ^ n order to avoid some bias 
in our results we follow Ref.[2J] and use the ratio CMB/BAO. This product cancels out some of the dependence 
on the sound horizon size at last scattering and a more model-independent constraint can be achieved in statistical 
analysis. In our analysis we use two ratio CMB/BAO measurements in redshifts z = 0.2 and 0.35, namely [24| (with 
one standard deviation error bars) 



d A {z 



D v (0.2) 
D V (Q.35) 



= 17.55 ±0.65, 
= 10.10 ±0.38 



(3) 

(4) 



where dx{z*) is the comoving angular-diameter distance to recombination, z* is the redshift of the last-scattering 
surface (z* = 1090) and Dy(z) is the so-called dilation scale. As emphasized in Ref. [24[ these two poins have a 
correlation coefficient of 0.337 14511 . In this point, we have to assume that brane models do not change the physics of 
the pre-recombination epochs [20[ . 



III. RESULTS 



A. Flat Universe 



The results of our statistical analyses to flat universe are shown in Figures la, lb and lc. Figure la shows contours 
of la and 2a on the H,\ — $7m plane when SDSS (MLCS2K2) and H(z) data points are considered separately (through 
this paper the priors 68 ± 2.8 and 73.8 ± 2.4 correspond to blue solid and red dash-dot curves, respectively). In this 
case Q\ can be find by using simply f^A = 1 — Q\ — &M- By considering two free parameters, we can find from the 
SNe la sample (la): fl M = 0.37 ± 0.11 and A = 0.0 + 0.022 ( X 2 = 237.84). On the other hand, from H(z) data 
points we find: D, M = 0.32 ± 0.1 and fi A = 0.0 + 0.006 {\ 2 = 15.87), and fl M = 0.24 ± 0.12 and Q x = 0.0 + 0.01 
(x 2 = 16.9) when the priors 68 ± 2.8 km/s/Mpc and 73.8 ± 2.4 km/s/Mpc are used on Ho, respectively. Note that no 
significant conflict between the results from H(z) priors is found and the bounds on £7 A from 21 H(z) data points are 
tighter than ones from 288 SNe la. Moreover, the contours from SNe la and H(z) analyses are almost ortogonal each 
other and a joint analysis can be used to impose tighter limits on the space parameters. 

Figure lb shows our results from a joint analysis. We find Ojvf = 0.37 ± 0.05 and f2 A = 0.0 + 0.0033 (\ 2 — 255.45) 
and n M = 0.33 ±0.04 and fl x = 0.0 + 0.0031 (\ 2 = 260.45) by using SNe la + H(z) with the priors 68 ±2.8 km/s/Mpc 
and 73.8 ± 2.4 km/s/Mpc on Hq, respectively. 
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FIG. 2: Confidence contours on the plane $1a — £Im to flat universe. Fig. 2a) Contours are drawn to la and 2a. The black 
dots corresponds to bounds from SNe la. Through this paper the H(z) priors 68 ± 2.8 and 73.8 ± 2.4 correspond to blue solid 
and red dash-dot curves, respectively. Fig. 2b) shows the contours to la and la from our joint analysis. Fig. 2c) displays the 
likelihoods to Q\. 

Panel lc display the likelihood for the H, x parameter, marginalizing on 51m- In this case we obtain tt x = 0.0 + 0.0018 
to both priors. 

B. Non-Flat Universe 

In this section, we allow deviations from flatness to prove the robustness of the previous Q x constraints using SNe 
la and H(z) data points. The results of our statistical analyses are shown in figures 2a, 2b and 2c (here we marginalize 
over Ha)- Figure 2a show contours of la and 2a on the Q, x — Q,m plane when SDSS (MLCS2K2) and H(z) data points 
are considered separately. In this case we can not constrain VLm and Q x simultaneously and, for instance: f2 m = 0.50 
and J1a = 0.005, and Ojvf = 0-20 and il x = 0.035 are permitted with high degree of confidence by using the SNe 
la sample. Furthermore, Q x = is excluded at 2a to Am = 0.1 and it is permitted at la to Qm = 0.3. Similar 
conclusion can be done to results from H(z). According with these comments a joint analysis is necessary in order to 
break the degenerescency on the (£l x , Qm) plane. 

Figure 2b displays the results by using a joint analysis involving SNe la and H(z) data points. We find (la): 
fl M = 0.38 ± 0.09, Q x = 0.0 + 0.0012 ( X 2 = 255.33) and Q M = 0.41 ± 0.08, fl x = 0.0 + 0.001 ( X 2 = 260.21) with the 
priors 68 ± 2.8 km/s/Mpc and 73.8 ± 2.4 km/s/Mpc on Hq, respectively. On the other hand, by marginalizing on il m 
we find: Q A = 0.65 ± 0.17 and Q A = 0.80 ± 0.18, respectively. 

Panel lc display the likelihood for the £l\ parameter, marginalizing on Qm and Qa- In this case we obtain (la) 
n x = 0.0 + 0.0018 and fi A = 0.0 + 0.0016 to 68 ± 2.8 km/s/Mpc and 73.8 ± 2.4 km/s/Mpc H priors, respectively. 
Therefore, the limits on Q\ are equivalent to the flat case and we see that the geometry has a negligible influence on 
estimates for this combination of data. 

Finally, in order to obtain tighter limits on space parameters we add in our analyses two points of the so-called 
ratio CMB/BAO. In figures 3a, 3b and 3c we show our results. We find (at la): Cl M = 0.28 ±0.03, fl x = 0.0 + 0.0001 
(X 2 = 262.88) and fl M = 0.28 ± 0.03, Q x = 0.0 + 0.00006 (x 2 = 264.21) with the priors 68 ± 2.8 and 73.8 ± 2.4, 
respectively. By marginalizing on $7m we find: ft a = 0.65 + 0.10 and 57a = 0.68 + 0.10 for each case, respectively. It is 
very important to stress that the limits derived by using SNe la + H(z) + CMB/BAO are more cosmological model 
independent than previous analyses where the A and R quantities were used [20(. The likelihood curves are plotted 
in figure 3c. We obtain fl x = 0.0 + 0.00004 and A = 0.0 + 0.00002 (la) to 68 ± 2.8 and 73.8 ± 2.4 priors, respectively. 



IV. CONCLUSIONS 



In this paper, we have studied observational bounds on cosmologies based on the flat and non-flat RSII brane 
models. Previous papers showed that SNe la data do not provide strong constraints on the parameters of RSII 
model. Moreover, the best fits did not give realistic values unless some priors to the matter and curvature densities 
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FIG. 3: Confidence contours on the plane Q\ — Qm to non-flat universe by using SNe la + H(z) + BAO/CMB ratio. Fig. 2a) 
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were assumed. However, as we have shown here, when we consider the joint analysis of 288 SNe la (SDSS by using 
MLCS2K2 method) and 21 H(z) measurements we obtain consistent results with other observations [47[ for matter 
density without imposing such priors. Moreover, the constraints on Cl\ are weakly dependent on the curvature 
parameter and Ho priors (68 ± 2.8 and 73.8 ± 2.4 km/s/Mpc), we have obtained at la: £l\ rs 0.0 + 0.002 to flat and 
non-flat Universes (see Figs, lc and 2c). It is important to emphasize that in the near future it will be possible to 
determine around 1000 H{z) values of the Hubble parameter at a 15% accuracy level [32, |48[. This fact shows that 
upcoming H(z) measurements plus SNe la may become a useful tool to impose constraints on Q\. 

We also have added the so-called CMB/BAO ratio in our analysis in order to obtain tighter limits on £1\. As 
is known, this ratio is more model independent than analyses done by using the A and R quantities which have 
parameters close to standard wCDM in their derivations. The joint analysis (SNe Ia+H(z)+CMB/BAO) gives the 
following values at la: fl x = 0.0 + 0.00004 and fl = 0.0 + 0.00002 with the priors 68 ± 2.8 and 73.8 ± 2.4, respectively. 
Finally, it is important to stress that the constraints for the brane tension found here are much weaker than the upper 
limit put by the analysis of Big Bang Nucleosynthesis (BBN) (Q\ < 10~ 20 ) [2(J, [4!|, however, our analysis is based on 
different physical phenomena and on physical processes that happened at a different cosmological epoch as compared 
to BBN and, therefore, constitute an independent test that must be performed in order to check the consistency of 
the model at large scale. 
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